Pressure-induced volume-collapsed tetragonal phase of CaFe2As2 
as seen via neutron scattering 
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Recent investigations of the superconducting iron-arsenide families have highhghted the role of 
pressure, be it chemical or mechanical, in fostering superconductivity. Here we report that CaFe2 As2 
undergoes a pressure-induced transition to a non-magnetic, volume "collapsed" tetragonal phase, 
which becomes superconducting at lower temperature. Spin-polarized total-energy calculations on 
the collapsed structure reveal that the magnetic Fe moment itself collapses, consistent with the 
absence of magnetic order in neutron diffraction. 
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Two recently discovered [H [3, 13 series of high tran- 
sition temperature (high- Tc) superconductors originate 
from the parent systems i?FeAsO {R — rare earth) and 
AFe2As2 {A — alkahne earth metal), which are tetrago- 
nal at room temperature but undergo an orthorhombic 
distortion in the range 100-220 K that is associated with 
the onset of antiferromagnetic order [H [gI. ItI. I^Js jJl^ fill . 
Tu ning the system via element substitution[2, y, Ij, Il2l . 
13, 14 1 or oxygen deficiency [isl. [l6j suppresses the mag- 



netic order and structural distortion in favor of super- 
conductivity (Tc's up to 55 K), with an overall behavior 
strikingly similar to the high-Tc copper oxide family of 
superconductors. 



The recent report 17| of pressure-induced supercon- 
ductivity in the parent CaFe2As2 compound opens an 
alternative path to superconductivity. Pressure sup- 
presses the distinct resistivity signature of the high- 
temperature structural and magnetic phase transition 
from 170 K at ambient pressure [3| to 128 K at 
0.35 GPa 17]. Superconductivity emerges with up to 
12 K for pressures between 0.23 GPa and 0.86 GPajnl- 
The pressure-induced superconductivity in CaFe2As2 
was confirmed 19] and followed by observations of su- 
perconductivity for BaFe2As2 and SrFe2As2 at signifi- 
cantly higher pressures [ioj. In CaFe2As2, a second high- 
temperature phase transition is observed above 0.55 GPa 
and 104 K by anomalies in the resistivity [17]. However, 
the nature of the phase at temperatures below this tran- 
sition and its relation to the ambient-pressure tetrago- 
nal, orthorhombic and pressure-induced superconducting 
phases are as yet unknown. 

Neutron scattering experiments on CaFe2As2 were 
performed to elucidate these issues. Special attention 



was paid to maintain experimental conditions closest to 
the reported macroscopic measurements and under well- 
defined hydrostatic pressure. Therefore, the experiments 
were performed on a polycrystalline sample prepared 
out of approximately 1.75 grams of single crystalline 
CaFe2As2 material grown using the procedure described 
in references and 21 1. The temperature profile for 



preparing this material was slightly modified (heating to 
1100°C and cooling over 50 hours to 600°C) to inhibit 
the formation of the reported[l^ needle-shaped impurity 
phase. Temperature-dependent resistance measurements 
on these crystals reproduced the data presented in refer- 
ences [l3| and [3] . The single crystals ( 300 pieces) were 
loaded with attempted random orientation into a He-gas 
pressure cell (maximum pressure 0.63 GPa) and mea- 
sured on the BTl high-resolution powder diffractometer 
at the NIST Center for Neutron Research. To reduce the 
effects of preferred orientation, the sample was oscillated 
over an angle of 36 deg during each measurement. For 
the temperature-dependent studies, the temperature was 
slowly changed with a maximum rate of 5 K/min while 
the pressure was adjusted and allowed to equilibrate be- 
tween measurements. 

Figure 1 shows neutron diffraction scans taken 
through the nuclear (0 2), (2 2 0)t, and magnetic 
(12 l)oR, magnetic diffractiou peaks at selected temper- 
atures and pressures. At 50 K and ambient pressure 
(A), the splitting of the (2 2 0)t into the orthorhombic 
(4 0)or/(0 4 0)oR peaks signals the transformation to 
the orthorhombic phase (Fig. 1(b)). This, together with 
the observation of the magnetic (1 2 1)or, magnetic peak 
(Fig. 1(c)), is consistent with previous x-ray and neutron 
diffraction measurements at ambient pressure [3, [S]- 
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FIG. 1: (color online) Scans through (a-b) nuclear and (c) 
magnetic peaks in neutron diffraction pattern at selected tem- 
peratures and pressures. Note that the diffraction peaks 
change position dramatically due to the significant changes in 
the lattice parameters. In (c), the magnetic (1 2 1)or, magnetic 
diffraction peak for point A is clearly observed above the back- 
ground taken at B. No new magnetic peaks for B and C were 
observed. Unlabelled peaks in the pattern arise from phases 
other than CaFe2As2, such as minor contamination from the 
Sn flux or Si02 (silica wool or pieces of the silica ampule from 
the single crystal growth), or the pressure cell. The subscripts 
denote the crystal structure used for indexing (T = tetrago- 
nal; OR = orthorhombic) . The offset between every data set 
is 200 Counts/3xl0* monitor in (a) and 300 Counts/SxlO* 
monitor in (b), respectively. 



Upon increasing pressure at T = 50 K, the structure 
remains orthorhombic and antiferromagnetic up to ap- 
proximately 0.24 GPa. 

Between 0.24 and 0.35 GPa, dramatic changes take 
place in the measured diffraction patterns. At pressures 
above 0.35 GPa (Fig. 1(c)), the magnetic peak is absent. 
No other magnetic magnetic peaks (e.g. corresponding to 
the AFl magnetic phase proposed in reference (i^l) are 
observed. The orthorhombic structure has transformed 
to a tetragonal phase, similar to the high-temperature 
ambient pressure structure, but with extraordinarily dif- 
ferent lattice parameters. This is most evident from the 
strong shift in the positions of the (0 2) and (2 2 0)t 
peaks at (B) in Figs. 1(a) and 1(b), respectively. The 
structure of this pressure-induced "collapsed" tetragonal 
phase is unchanged in the superconducting state deter- 
mined by measurements at 4 K and under 0.48 GPa. 

The central region (shown in yellow) of Figs. 2(b) 
and (c) shows the results of Rietveld refinements of 
the lattice parameters and volume for the "collapsed" 
tetragonal phase. The structure data (lattice parameters 
and atomic-position parameter zas of As) were deter- 
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FIG. 2: (color online) Schematic phase diagram and pressure 
and temperature dependence of the lattice parameters, unit 
cell volume and As-Fe-As bond angles, (a) (left panel) Lines 
in the schematic p — T diagram denote the high-temperature 
and superconductivity phase lines determined in reference 
[13]. Points A-D label the pressures and temperatures for 
the diffraction data shown in Figure 1. (a) (right panel) The 
unit cell of the tetragonal phase of CaFe2As2. (b-d) Pressure 
dependence (left panels) at T = 50 K and temperature de- 
pendence (right panels) at p = 0.63 GPa (b) of the lattice 
parameters, (c) of the unit cell volume and c/a ratio, and (d) 
of the As-Fe-As bond angle. 



mined by Rietveld refinements using the GSAS software 
package [23j. We find an astonishing 9.5% reduction in 
the c-lattice parameter with respect to the orthorhombic 
phase and a nearly 5% decrease in the unit cell volume. 
Even more striking is the reduction of the c/a ratio, a 
key parameter for bond geometries in the iron arsenides, 
by nearly 11%. As a consequence, the As-Fe-As bond 
angles change strongly as illustrated in Fig. 2(d). 

With the pressure maintained at 0.63 GPa, the temper- 
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ature was raised in 50 K steps (right panels of Fig. 2) . Be- 
tween 150 and 200 K an isostructural transition between 
the low-temperature "collapsed" tetragonal phase and 
the high-temperature tetragonal structure is observed. 
Upon release of the pressure at 250 K, the curves labelled 
(D) in Fig. 1 show only small changes in the lattice pa- 
rameters between 0.63 GPa and ambient pressure, pro- 
viding a measure of the modest, but strongly anisotropic 
compressibility of the high-temperature phase. 

We note that there is a difference of about 50 K be- 
tween the temperature of the isostructural transition at 
0.63 GPa measured here and that reported in trans por t 
measurements [l7||. However, as pointed out in Ref. 
the resistive anomalies are rather broad in applied pres- 
sure, and different criteria for the definition of transition 
temperatures can shift temperature assignments. In ad- 



dition, the data in Ref. 17| were taken with decreas 



ing temperature whereas here the temperature was step- 
wise increased. With these uncertainties understood, the 
tetragonal-to- "collapsed" tetragonal transition appears 
to be responsible for the loss of resistivity whose locus de- 
fines the high-temperature high-pressure phase line found 
in Ref. (ITJ and shown in Fig. 2(a). 

In order to relate the volume change to relative changes 
in the unit cell dimensions, and to verify the stability of 
this phase, spin-polarized total-energy calculations were 
performed for volume changes of AV/V = 0% (for am- 
bient pressure) and AV/V — —5% (for the "collapsed" 
phase). The local density approximation was employed, 
using the full potential linearized augmented plane wave 
method with the Perdew-Wang 1992 functionalfi^]. The 
precision of the total energy is 0.01 mRyd/cell, much 
smaller than the size of the symbols in Fig. 3. 

From the blue (dark grey) curves in Fig. 3(a) we see 
that, for ambient pressure, the orthorhombic magnetic 
phase is lowest in energy, consistent with our ambient- 
pressure low-temperature measurement. The red (light 
grey) curves in Fig. 3(a) show that the tetragonal phase is 
lowest in total energy for the 5% volume reduction. The 
minimum energy of this "collapsed" tetragonal phase is 
found at c/a ~ 2.65, close to the experimental value of 
2.67 (Fig. 2(c)). 

The c/a-dependence of the spin-polarized total-energy 
calculations for the "collapsed" phase can be corre- 
lated with a loss of the Fe magnetic moment, as shown 
in Fig. 3(b). Both the spin-polarized and non-spin- 
polarized calculations yield the same total-energy min- 
imum for the non-magnetic "collapsed" tetragonal phase 
at the same c/a ratio. The astonishing result of a 
quenched magnetic moment ground state is consistent 
with our experimental observation of the loss of mag- 
netic order in the "collapsed" tetragonal phase. The band 
structure calculations also indicate that several bands 
cross the Fermi level at the pressure-induced transition. 

The principal result of these neutron diffraction mea- 
surements is the discovery of a transition from the 
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FIG. 3; (color online) Summary of the results of the neu- 
tron diffraction measurements and total-energy calculations, 
(a) Spin-polarized (SP) and non-spin-polarized (NSP) total- 
energy calculations for AV/V = 0% and AV/V = -5% for 
the tetragonal (T) and orthorhombic (OR) phases (Note the 
different energy scales). The tetragonal notation is used for 
the c/a ratio, (b) For the "collapsed" phase, the Fe moment is 
quenched at the minimum in total energy in the spin-polarized 
calculation, (c) Schematic p — T diagram based upon the 
diffraction and transport[l3] measurements. Dashed lines rep- 
resent estimates of the phase boundaries. The cross hatched 
area indicates a region where precise details of how these tran- 
sition lines intersect still needs to be determined, (d) Corre- 
lation between the superconducting transition temperatures 
and the variance of the As-Fe-As bond angle. The data 
were obtained from Refs. f5| and [H, |M 113, ill, El] for the 
iiFeAs(0/F) compounds and from Ref [J] for (Ba/K)Fe2As2, 
respectively. The data for CaFe2As2 are derived from the 
present measurements. 



magnetically ordered orthorhombic phase to a non- 
magnetically ordered "collapsed" tetragonal phase pre- 
ceding the onset of superconductivity. Further, the 
second, higher pressure, transition noted in transport 
measurements |17| has been identified as an isostructural 
transition between the pressure-induced "collapsed" 
phase and the high-temperature tetragonal structure. 
The observed volume reduction can, for example, serve to 
increase the charge-carrier density. The schematic phase 
diagram in Fig. 3(c) summarizes our findings. Our results 
show that the pressure- induced superconductivity [l3| 
in CaFe2As2 emerges from the "collapsed" tetragonal 
phase rather than the magnetically ordered, orthorhom- 
bic phase or the high-temperature tetragonal phase. 

Anomalous changes in the unit cell volume and lattice 
constants have also been noted in the i?FeAsO system. 
In superconducting fluorine- free oxygen-deficient samples 
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of NdFeAsOi_5, a surprising discontinuous decrease in 
the lattice parameters and unit cell volume {/S.V/V — 
— 1.8%) was found for 5 = 0.4, where a maximum in the 
superconducting volume fraction is observed Fur- 
thermore, a correlation between superconducting transi- 
tion temperatures for i?FeAsOi_5 and the unit cell di- 
mensions is reported[l^. In light of our results, it ap- 
pears that chemical substitution and the introduction of 
oxygen deficiency likely play a dual role in the iron ar- 
senide superconductors by increasing the charge carrier 
density and changing the "chemical pressure" . It is not 
yet clear which has the greater impact upon supercon- 
ductivity. 

Given that the i?FeAsO and ^Fe2As2 families share a 
common structural element (FeAs layers) and similar pre- 
requisites for superconductivity (e.g. suppression of mag- 
netic order) it is useful to elucidate structural quantities 
that are shared and can be correlated with superconduct- 
ing properties. Of particular interest is the comparison of 
the As-Fe-As bond angles (25| with the ideal tetrahedral 
value of 109.47 deg. For this special value all Fe atoms 
are coordinated in ideally formed tetrahedrons with iden- 
tical high-symmetric Fe-As-Fe bonding geometries [high- 
lighted in brown in the right panel of Fig. 2(a)]. This 
symmetry is broken for As-Fe-As bond angles deviating 
from this ideal angle yielding two or three different values 
in the tetragonal or orthorhombic structure, respectively. 
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The variance cr^ — ^X)[(t'ond angle) — 109.47 deg]^ in 
1 

deviation of the Fe-As-Fe bonding angles from the ideal 
value parameterizes the strength of this symmetry break- 
ing. 

In Figure 3(d) we plot the measured superconduct- 
ing transition temperatures Tc, for those iron-arsenide 
compounds that are superconducting, as a function of 
this variance cr^ of the As-Fe-As bond angles. As in- 
creases, corresponding to greater deviations in the bond 
angle from the ideal tetrahedral angle, Tc decreases. The 
pressure-dependent properties of the "collapsed" tetrag- 
onal phase of CaFe2As2 described above clearly continue 
this trend and allow generalizing our result to the family 
of FeAs based superconductors. The high value for in 
the "collapsed" tetragonal phase for CaFe2As2 is consis- 
tent with the low Tc in comparison to other FeAs based 
superconductors. The observed correlation between Tc 
and cr^ points to the importance of the structure and 
the symmetry in the FeAs network for the superconduct- 
ing state. Properties that are sensitive to the As-Fe-As 
bonding geometry and its symmetry, such as anisotropic 
magnetic or elastic couplings in the FeAs network, seem 
strongly involved in the superconducting pairing. 

These results highlight intriguing questions that point 
to the potential complexity of the superconducting state 
in the iron arsenides. From the analysis of Tc as a func- 
tion of the As-Fe-As bond angle variance, we have found 
that superconductivity in the "collapsed" tetragonal 



phase fits well within the general trend observed for 
the doped iron arsenides, implying a common supercon- 
ducting pairing mechanism. It has been suggested that 
spin fluctuations are responsible for the electron pairing 
in this class of superconductors 0. The apparent loss 
of a static moment in the "collapsed" tetragonal phase 
may seem inconsistent with such a magnetic pairing 
mechanism for the pressure-induced superconductivity 
in CaFe2As2. However, the strong pair-breaking effect 
of local moments is eliminated and superconductivity 
mediated by paramagnons[3lj remains a possibility. 
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Basic Energy Sciences under Contract No. DE-AC02- 
07CH11358. M. S. Torikachvili acknowledges support 
by the National Science Foundation under DMR-0306165 
and DMR-0805335. 



APPENDIX A: DETAILS OF THE NEUTRON 
DIFFRACTION MEASUREMENTS AND THEIR 
ANALYSIS INCLUDING STRUCTURE DATA 

The neutron diffraction measurements were performed 
on the high resolution powder diffractometer BTl at the 
NIST Center for Neutron Research using a wavelength of 
2.0782 A selected by a Ge (3 1 1) monochromator. The 
collimation of the incident beam was set at 15'. 

Special attention was paid to maintaining experimen- 
tal conditions closest to the reported macroscopic mea- 
surements and to perform the study under well-defined 
hydrostatic pressure. The soft and ductile nature of 
the compound presents challenges for powder diffraction 
measurements since flux-grown single crystals tend to 
smear and shear when ground into a powder. The effects 
of grinding are clearly observed as broadened peaks in x- 
ray powder diffraction measurements, and the modifica- 
tion of physical properties associated with grinding have 
not yet been characterized. Therefore, as grown single 
crystals ('^500-1000) were loaded into an Al-alloy He-gas 
pressure cell and cooled using a top-loading closed-cycle 
cryogenic refrigerator. The pressure cell was connected 
to a pressurizing intensifier through a high pressure cap- 
illary. Hydrostatic pressure was maintained throughout 
the measurements since the temperature was kept well 
above the melting curve for helium (for pressures up to 
0.63 GPa). For the temperature dependence studies, the 
temperature was slowly changed with a maximum rate 
of 5 K/min while the pressure was adjusted and allowed 
to equilibrate between measurements. 

To reduce the effects of preferred orientation, the sam- 
ple was oscillated over an angle of 36 deg during each 
measurement. Nevertheless, a degree of preferred ori- 
entation remained but was adequately modelled in sub- 
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FIG. 4: (color online) Rietveld analysis of the neutron diffrac- 
tion pattern at p = Gpa and T = 50 K. Areas with strong 
contributions from the pressure cell are excluded. 
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FIG. 6: (color online) Rietveld analysis of the neutron diffrac- 
tion pattern at p = 0.63 Gpa and T = 250 K. Areas with 
strong contributions from the pressure cell are excluded. 
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FIG. 5: (color online) Rietveld analysis of the neutron diffrac- 
tion pattern at p = 0.63 Gpa and T = 50 K. Areas with strong 
contributions from the pressure cell are excluded. 



sequent Rietveld refinements using the GSAS software 
package. For each diffraction pattern approximately 70 
reflections were used to refine the lattice parameters, the 
z coordinate of the As ions and the 12 parameters associ- 
ated with corrections for preferred orientation. We point 
out that the correction for preferred orientation remained 
constant over all pressures and temperatures measured 
and were taken as constants in the fits. Typical R- values 
(wRp), representing the goodness-of-fit, were between 4- 
5% for all fits demonstrating the accuracy of the model 
employed. 



APPENDIX B: DETAILS OF THE 
SPIN-POLARIZED TOTAL ENERGY 
CALCULATIONS 

For the spin-polarized and non-spin-polarized cal- 
culations, the local density approximation was em- 
ployed, using the full potential linearized augmented 
plane wave method with the Perdew-Wang 1992 func- 
tional. The convergence criterion for the total energy 
was 0.01 mRyd/cell. The calculations were done in the 
tetragonal phase(two formula units) for two different cell 
volumes as the c/a ratio was varied: 0% volume reduc- 
tion using the experimentally determined lattice param- 
eters and for a 5% volume reduction. The calculations 
were also performed for the orthorhombic phase, with 
ajh = 1.02, for both 0% and 5% volume reduction. The 
experimental value of the internal parameter was used for 
the calculations. The i?MT*-f^^max that determines matrix 
size (the number of the basis functions), where i?MT is the 
smallest of all atomic sphere radii and if^max is the plane 
wave cut-off was set to 8.0. The muffin-tin radii (Rwr) 
for the 0% volume reduction calculation were 2.2, 2.1 and 
2.1 atomic unit for Ca, Fe and As, respectively. For the 
5% reduction calculation, the i?MT also were scaled to 
keep the matrix size the same. There were 100 fc points 
used in the irreducible Brillouin zone. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 
J. Am. Chem. Soc. 130, 3296 (2008). 

[2] H. Takahashi, K. Igawa, K. Arii, Y. Kamihara, M. Hi- 
rano, and H. Hosono, Nature (London) 453, 376 (2008). 

[3] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and 
D. F. Fang, Nature (London) 453, 761 (2008). 

[4] M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 
101, 107006 (2008). 



6 



TABLE I: Results of Rietveld refinements for pressure-dependent measurements at a temperature T = 50 K 



Pressure p (GPa) 





0.115 


0.24 


0.35 


0.47 


0.63 


a(A) 
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h{X) 
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5.4577(2) 


5.4575(2) 
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1 A ar>'-70k(ci\ 
lU.Do lijyO) 
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iU.DU i o\i ) 


Volume V(A^) 


352.68(4) 


351.94(3) 


351.11(3) 


168.44(1) 


168.07(1) 


167.85(1) 


ZAs 
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Fe-Fe (A) (2x) 


2.7288(1) 


2.72882(7) 


2.72874(8) 








Fe-Fe (A) (4x) 








2.8137(9) 


2.8132(1) 


2.8128(1) 



TABLE IL Results of Rietveld refinements for pressure-dependent measurements at a pressure p = 0.63 GPa 



Temperature T (K) 


4 (p = 0.47 GPa) 


50 


100 


150 


200 


250 


250 (p = GPa) 


a{A) 


3.9797(3) 


3.9780(1) 


3.9760(3) 


3.9724(2) 
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3.8944(2) 


3.8915(1) 


c(A) 
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10.6073(7) 


10.633(1) 


10.683(1) 


11.438(2) 


11.530(2) 


11.690(1) 


Volume V(A^) 


168.33(3) 


167.85(1) 


167.95(3) 


168.57(3) 


174.11(3) 


174.87(3) 


177.04(2) 


ZAs 


0.365(2) 


0.3663(5) 


0.3668(9) 


0.3669(9) 


0.367(1) 


0.3652(8) 


0.372(1) 


Fe-As (A) 


2.333(9) 


2.340(3) 


2.344(5) 
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2.365(9) 


2.357(5) 


2.410(9) 


Fe-Fe (A) (4x) 
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2.8128(1) 
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2.8089(2) 


2.7588(2) 


2.7538(1) 


2.7517(1) 
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